home *** CD-ROM | disk | FTP | other *** search
/ The Atari Compendium / The Atari Compendium (Toad Computers) (1994).iso / files / prgtools / programm.ing / falcon / nt_dsp1.lzh / NT_DSP1.MSA / FNTNS / SQRT1.ASM < prev    next >
Encoding:
Assembly Source File  |  1989-01-24  |  1020 b   |  30 lines
  1. ;
  2. ; This program originally available on the Motorola DSP bulletin board.
  3. ; It is provided under a DISCLAMER OF WARRANTY available from
  4. ; Motorola DSP Operation, 6501 Wm. Cannon Drive W., Austin, Tx., 78735.
  5. ; 
  6. ; Square Root by Polynomial Approximation, 7-bit Accuracy.
  7. ; 
  8. ; Last Update 26 Jan 87   Version 1.0
  9. ;
  10. sqrt1   macro
  11. sqrt1   ident   1,0
  12. ;
  13. ;       Does a square root by polynomial aproximation, 7 bit accuracy.
  14. ;       sqrt(x)=-.2831102*x*x+2*.5031642*x+.272661
  15. ;                  a2             a1          a0
  16. ;       valid for:  .25<= x < 1.0
  17. ;
  18. ;       input value in x0, output in regester A.
  19. ;
  20. ;       r1 initially points to the coefficients in y memory in
  21. ;       this order: a1, a2, a0.
  22. ;
  23. ;
  24.         mpyr    x0,x0,a  y:(r1)+,y0 ;x**2, get a1
  25.         mpy     x0,y0,a  a,x1           ;x*a1, save x**2
  26.         mac     x0,y0,a  y:(r1)+,y0     ;2*x*a1, get a2
  27.         macr    x1,y0,a  y:(r1)+,y0     ;a2* x**2, get a0
  28.         add     y0,a                    ;add in a0
  29.         endm
  30.